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Abstract 

A novel tracking paradigm for flying geometric trajectories using teth- 
ered kites is presented. It is shown how the differential-geometric notion 
of turning angle can be used as a one-dimensional representation of the 
kite trajectory, and how this leads to a single-input single-output (SISO) 
tracking problem. Based on this principle a Lyapunov-based nonlinear 
adaptive controller is developed that only needs control derivatives of the 
kite aerodynamic model. The resulting controller is validated using sim- 
ulations with a point-mass kite model. 

1 Introduction 

In 1980, Loyd wrote a seminal paper exploring the possibility of generating 
electrical power using the pulling force of tethered airfoils, i.e., kites [19]. During 
the oil glut of the 1980s, however, interest in wind energy dropped [1]. It was 
only during the turn of the 21 st century that interest in kite power (and traction) 
picked up again with the work of Meijaard, Ockels, Schwab and Diehl [20, 21, 8]. 

Currently, almost the whole wind energy market makes use of horizontal 
axis wind turbines. Wind turbines harvest the energy relatively close to the 
ground. But realizing that up to a height of approximately 10 kilometres the 
wind velocity increases with altitude [18], and that the wind power raises with 
the wind speed to the third power [19], one concludes that there is a large 
potential to extract wind energy at higher altitudes. A research group at the 
Institute for Applied Sustainable Science, Engineering and Technology (ASSET) 
at Delft University of Technology is investigating a wind energy converter using 
kites as lifting surfaces, i.e., the Laddermill [21]. 

The design of automatic controllers for kite power systems has been a subject 
of ongoing research [25, 26, 15, 12, 13, 9, 7]. In this paper we discuss two 
contributions. First, we illustrate a trajectory tracking concept based on the 
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angle between the velocity vector and the horizontal. This concept is restricted 
to the motion of the kite in directions perpendicular to the vector connecting 
the earth tether attachment point and the kite, leaving control of the reeling 
in and reeling out of the tether up to a separate controller. To get a first 
intuitive picture of this idea, let us for a moment project the crosswind figure- 
eight trajectory of a kite onto the plane perpendicular to the wind vector. As a 
kite traces out a figure eight on this plane, the angle the velocity vector makes 
with the horizontal oscillates between certain minimum and maximum values; 
see Figures 1 and 2. 




Figure 1: Tangent vectors along path. Figure 2: Tangent vectors translated 

to the origin. 

Instead of measuring this angle directly we use the more general notion of 
turning angle [11] from elementary differential geometry (see the Appendices). 
This angle is not restricted to [— tt, tt] and it is continuous. It is cumulative in 
the sense that multiple turns around the origin of the velocity vector will result 
in a multiple of 2ir, but winding backwards will reduce the turning angle by the 
same amount for every turn. This property makes the turning angle well suited 
for control purposes. By considering the arc length as independent variable, the 
turning angle - together with initial conditions - becomes a representation of 
the image of a curve that is independent of parametrization. 

Secondly, we derive an adaptive turning angle tracking controller that does 
not need a full aerodynamic kite model. Instead it only uses the control deriva- 
tives relating control input increments to steering force increments; the rest of 
the required data is measured directly. We provide a stability proof for the 
resulting controller and demonstrate its performance with simulations. 

The remainder of this paper is organized as follows. In the next section we 
discuss the physics of the kite control problem and touch on previous research. In 
section 3 we formalize a geometric model of the dynamics of the kite trajectory, 
and in section 4 we formulate and solve the control problem. We illustrate 
the performance of the control law with simulations in section 5. Finally, the 
appendices list relevant results from the field of differential geometry. 
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2 Motivation 

The Laddermill uses kites to pull a tether from a drum driving a generator 
to produce electricity. The operational mode consists of two distinct phases: 
a reeling-out phase pulling the tether from the drum producing power and a 
reeling-in phase which consumes energy; see Figure 3. It is desired to control 
the kite in such a manner that a net power production remains between the 
reeling-out and reeling-in phase. By flying the kite in directions perpendicular 
to the wind, the magnitude of the airflow along the airfoil is increased [19]. This 
principle is comparable to the crosswind rotation of wind turbine blades. 




Figure 3: The Laddermill during the reel-out phase [7]. 

The control mechanism developed at the ASSET institute consists of a con- 
trol pod suspended under the kite. The control pod is connected to the drum on 
the ground using a single tether and four bridle lines lead from the control pod 
to the kite. The pod contains two servo motors, one controlling the length dif- 
ference of the two back lines, the steering lines, and the other motor controlling 
the overall length of the steering lines. By changing the length difference of the 
steering lines the kite deforms; see Figure 4. The aerodynamics in turn generate 
steering yaw moments. A detailed explanation is available in [5]. Regulation of 
the total length of the steering lines provides some control over the attack angle 
of the kite itself. This is known as the powering and the depowering of the kite. 




Figure 4: Kite before and after applying a steering input [5]. 

Previous research into automatic control of kites has focused on trajectory 
tracking using linearization [25, 26] and Model Predictive Control (MPC) [15, 
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12, 13, 9]. All of these approaches require a fully validated kite model that can 
be applied in real time. Multi-body models such as those developed in [5] are 
precise but too computationally intensive for control design, and condensation 
of kite dynamics into a rigid body model using system identification techniques 
has proven to be difficult [7]. 

3 Geometric modeling and system 

We are interested in controlling flight in directions perpendicular to that of the 
radial vector r connecting the earth tether attachment point and the kite. The 
control of the tether and of the power lines is left to a separate controller. In 
order to derive a control system that is invariant of the tether length, we project 
the flight trajectory onto the upper half of the unit sphere centered on the earth 
tether attachment point. This yields the projected path 7 := r/||r||. 

3.1 Geometric model 

As mentioned in the introduction, we will control the turning angle of the kite 
trajectory. In order to be able to compute the turning angle on the upper half of 
the unit sphere we need a coordinate system. For this, note that the coordinate 
patch 

p : [0,27r) x [0,7r/2] — > R 3 , (v,w) h-> (cos w cos v, cos w sin v, sin w) . 

yields a continuous one-to-one (except for the pole) correspondence between the 
upper half of the sphere and the set [0, 2ir) x [0, n/2]. The vector fields dp/dv 
and dp/dw are orthogonal, and by normalizing them we obtain an orthonormal 
basis of the tangent plane, 

ei := (— sinu, cosw, 0), e 2 := (— cos w sin w, — sin u sin w, cosw), (1) 

as illustrated in Figure 12. 




Figure 5: Coordinate system on the upper half of the unit sphere. 
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The vector field ei has no vertical component, which intuitively makes it a 
suitable reference for the definition of a turning angle. To make this precise, we 
apply Liouville's theorem [11] for which we need the curvatures of the coordinate 
system. The coordinate curves w i-> p(v,w) are parts of great circles, whence 
their geodesic curvature (n g ) 2 is zero. The curves v i-> p(v,w), however, have 
geodesic curvature (/%)i = sinw. By Liouville's theorem, the turning angle of 
the projected kite trajectory with respect to the vector field ei is given by the 
integral equation 

9[i\ :=0[7,ei] = [ K g [>y] - (k„)i cos 0[ 7 , ei] ds + 0[ 7 , ei] , (2) 

where K g [y] is the geodesic curvature of the projected kite trajectory and 8[j, e^o 
is the initial value of the turning angle. The geodesic curvature of the projected 
kite trajectory can be calculated as 

xy—'xy 

with x and y the coordinates of the velocity of the projected trajectory, expressed 
in the (ei, e 2 )-basis of the tangent plane. 



3.2 Path length 

The projected velocity of the kite may sometimes have a component perpen- 
dicular to the target trajectory 7 t . If we zig-zag around the target trajectory 
and straighten out the traversed path we obtain a longer arc length than if we 
would straighten out the target trajectory. We need to take this difference into 
account, for otherwise we would move through the target trajectory quicker than 
intended. For this reason we introduce a corrected path length that disregards 
the perpendicular component. 

Let T t := cos6>[7 t ]ei + sin#[7 t ]e2, where (ei,e2) is the orthonormal basis 
of the tangent plane at 7 t given by (1). This unit vector field is tangent to 
the target trajectory. We transport this tangent vector to the tangent space 
at 7 by means of the operator R, which rotates 7 to 7 t along their connecting 
geodesic. Taking the inner product 7 • (R T T t ), we obtain the component of 
the projected velocity tangential to the target trajectory. This motivates the 
following definition. 

Definition 3.1. We refer to the integral 

s c := Jj • (R T T t ) dt 

as the corrected path length. 



3.3 Control system 

Our objective is to track the given target trajectory 7 t which, as hinted at 
previously, we will do by controlling the turning angle #[7]. We show that we 
can achieve this objective using a single steering control input. This way we 
construct a one-dimensional single- input, single-output tracking problem. 
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In order to arrive at a set of equations that are affine in their control input 
we base our discussion on steering input increments u := u — u* . Central to 
this linearizing approach is the following assumption. 

Assumption 3.1. The acceleration 7 is locally affine in feasible 1 control in- 
crements u. That is, 



<9 3 7 
dudt 2 



> 



<9 4 7 . 

11 

du 2 dt 2 



Linearizing the derivative of the turning angle in the control input u we find, 
using (2), (3), and the definition of the corrected path length, that 



d9[j] _ d9[~f] 



dt 



d 2 e\ 1 ] 



dt 

d9[-y] ds 



ds dt 
xy ~ 'xy 



+ 



dudt 
d 2 6[j] ds 



duds dt 



(4) 



-j-e z cos6»[7]||7|[ 



' du 



a^- 



where 7 • e z denotes the (earth) vertical component of the projected trajectory. 
In deriving our control loop in the next section we will make use of the direct 
coupling between the control increment u and the derivative of the turning 
angle. This principle is similar to the incremental nonlinear dynamic inversion 
approach for control of aircraft as developed by Sieberling ct al. [23]. 

For the derivatives dx/du and dy/du, representing the sensitivities of the 
accelerations to changes in the steering input, we need a model. Let r = ||r|| 
and note that 

.. d 2 



7 



rr 



2^+2 — . 



dt 2 \r ) r r^ r^ r° 
Expressed in the tangential frame the first two coordinates of the radial vector 
r are zero. As a result, and due to the fact that only the control derivatives 
of the accelerations are non-zero, we see that the sensitivities, expressed in the 
tangential frame, are given by 



dx pvlStC 
du 2mr 



■t 

x.u 



du 



2mr 



(5) 



Here C l x u and C* u are the unitless control derivatives expressed in the tan- 
gential frame, p is the density of air, v a is the airspeed, S t is the surface area 
of the tips and m is the combined mass of the kite and the control pod. Let 
C u — [C x . u Cy. u C ZtU ] denote the control derivatives expressed in a body-fixed 
reference frame. 

Kites deform during flight and at any point in time the deformation depends 
on the form the kite had previously. Hence the deformation of the kite is path- 
dependent. In order to model the control derivatives, we assume path-dependent 
deformation to have a - for control purposes - negligible influence on the steering 
aerodynamics. 



x We call a control input increment feasible when the actuator can implement it within a 
single sampling interval. 
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Assumption 3.2. The unitless control derivatives C u , expressed in the body- 
fixed reference frame, are a function of the current steering input u. 

Note that this assumption can be relaxed by including the dependence on 
additional state variables. 

We approximate the function u i-> C„, mapping states to control derivatives, 
using B-spline networks as in [24]. Let Wj denote the weights and b the basis 
functions of the B-spline network for the ith control derivative with i e {x, y, z}, 
so that the zth B-spline network is given by C\. u — w, • b. The number of basis 
functions used for the spline networks must be chosen in such a way as to avoid 
problems with under- or overparametrization. This can be done using a trial 
and error procedure. 

By continuously measuring the position, linear velocity, linear acceleration, 
airspeed and attitude of the kite, we obtain a nonlinear system of the form 

d M = mi ],t)+B(t)u. (6) 

The measurements are implicit in the time-dependence of the functions / and 
B, which are given by 



f = d9h] 



dt 
and 



*«_ 7 .e z cos 0[ 7 ]||7ll 



B 



dudt 



,dy_ _ dx- 



II7II 2 



By choosing our system in this way the need for a full aerodynamic kite model 
vanishes. 

From (4) and (5) we see that the function B is linear in the control derivatives 
C»,« = Wj ■ b, that is, 

B = J2 A * w * ' b ' ( 8 ) 

i 

for some coefficients Aj G R. This shows that the system (6) is affine in both the 
control input u as well as in the control derivative spline weights Wj , a property 
that will allow us to design an adaptive control loop with relative ease in the 
next section. 



4 Control 

In this section we design our control loop and control derivative estimators by 
means of a control-Lyapunov approach. Our approach is closely related to those 
described in [24, 3, 16, 17, 14]. 

Our objective is to track a given reference projected kite trajectory. To this 
end, we construct a flight controller consisting of an inner and an outer loop. 
The inner loop controls the turning angle of the kite trajectory 9[f] towards a 
target turning angle 9 t . The outer loop computes a target turning angle 9 t that 
steers the kite onto the target trajectory. 
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4.1 Inner loop 



Our inner loop cost function is the difference between the actual and desired 
turning angle: 

e:=0ky]-6 t . (9) 

In order to track the reference trajectory we must find a control law that renders 
e -> 0. 

The steering dynamics of a kite may subtly change over time due to weather 
patterns and wear. In order to automatically adapt to these changes we aim 
for an adaptive control approach that maintains and updates an internal rep- 
resentation of the steering control derivatives. We will denote this internal 
representation as Wj with i <G {x, y, z}. 

Summarizing, we formulate the following control problem: 

Problem 4.1. Find a steering control law and a control derivative B-spline 
weight estimator update laws that such that e — > as t — > oo. 



4.1.1 Control design 



In this section we propose steering control and estimator update laws to solve 
Problem 4.1 using a control-Lyapunov approach. 



Proposition 4.1. Consider the system (6). The control law 



deb] 

dt 



dth 
dt 



Ke 



(10) 



with gain K € R + , and the control derivative B-spline weight estimator update 
laws 

^ := reAifib, (11) 
with adaptation gain T € R + , render the Lyapunov function 



V := V , 

2 2 



(12) 



decreasing along the flow. 

Proof. The time-derivative of the Lyapunov function V is given by 



dV_ 
~dt 



de 



dwj 
~d~t 



(13) 



with, using (4), 



de _ d6[y] 
dt ~ dt 



d 2 e[ 1 ] 



dudt 



dt 



The estimator update laws (11) cancel out all terms with Wj. Indeed, replacing 
these into (13), we obtain, using (7) and (8), that 



dV 


\d9[~f] 




d0 t ~ 


— — = e 
dt 


dt 




dt 



+ euj^ A^Wi • b. 



(14) 



8 



Finally, substituting the control law (10) into (14), we see that 

dt 

The gain K € R + renders V strictly decreasing whenever e ^ 0. 



□ 



In the previous analysis we neglected rate and magnitude constraints of 
the control signal. Let u a denote the control increment implemented by the 
actuator. Since in general u a ^ u, the Lyapunov function (12) is no longer 
guaranteed to be decreasing and the parameter estimation process may diverge. 
In order to remedy this situation we introduce the modified tracking error as in 
[10], 

with 

d( ^ , d 2 9[-y] 
dt 



-KC 



dudt 



(u a - u), 



(15) 



providing a correction for the influence of the actuator limitations, where the 
second derivative is evaluated using the estimated weights w, . Note that in the 
absence of actuator constraints e m — > e. 

In order to render the modified tracking error to zero, we adjust the es- 
timator update laws (11) to be proportional to e m . The associated modified 
Lyapunov function is - by construction - decreasing along the flow, indepen- 
dently of whether or not u a = u: 

Proposition 4.2. Consider the system (6). Any control law and the modified 
control derivative B-spline weight estimator update laws 

dt 

render the modified Lyapunov function 



Te m \iU a b 



(16) 



V - l -e 2 



(17) 



decreasing along the flow. 

Proof. The time-derivative of the modified Lyapunov function (17) is given by 



dVrn 

dt 



dp 

'~d~t 



+ F- 1 ^(w,-w J ) 



~df 



with, using (15) and (4), 
de m „, , d 2 9[ 1 ] 



dt 



») 



dudt 



(u a -u) + 



Mb] 



U* ,Wi 



dt 



+ 



d 2 9[j} 



dudt 



(18) 



dfft 
dt ' 



The modified estimator update laws (16) cancel out all terms with w,. Indeed, 
replacing these into (18) we obtain, using (7) and (8), that 



dV m 

dt 



K(e 







d 2 9[ 7 ] 



dudt 



\d0h] 




d0 t ~ 


dt 


U* 


dt 



(u a - u) 

«*,w, 

+ e m u a ^2 ^ 



(19) 
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Noting that u a = (u a — u) + u and substituting the canonical control law (10) 
into (19) we see that 



dt 



= -Kei 



The gain K e R+ renders V m strictly decreasing whenever e m 7^ 0. 
4.1.2 Control performance 



□ 



In this section we investigate the performance of the proposed control and es- 
timator update laws. In the case that u a = u, the function V defined in (12) 
is decreasing along the flow. This implies convergence of the kite trajectory 
turning angle 9 [7] to the target turning angle 9 t ■ 

Theorem 4.1. Consider the system (6), the control law (10) and the estimator 
update laws (16). Assume that there is a t > such that for all t > t 

Il7(*)ll>0 

holds, so that for t > t the kite trajectory is a regular curve. Additionally 
assume that for all t > to the geometric condition 



dj 
du 



IK*), 



i.e. that the given two vectors are not parallel, is satisfied, so that the control 
law is well-defined. Finally, assume that for some M G R + and for all t > 
the bounds 



<M, \6[ 7 )(t) - 9 t (t)\ <M 



dOb] 




< M, 


dB t 




< M 


dt 


t 




dt 


t 





hold. Then, for any initial condition and for any reference trajectory the turning 
angle of the kite trajectory converges to the target: 

\im9[ 7 ](t) = e t (t). 

t— ^ 00 

Proof. Barbalat's Lemma [16] applied to the Lyapunov function V shows that 
e -)• 0. □ 

We proceed to show that, even when u a ^ u, the modified parameter esti- 
mation process is stable in the sense that the weight estimation errors remain 
bounded: 

Theorem 4.2. Consider the system (6) and the estimator update laws (16). 
Assume that there is a t > such that for all t > t the conditions 



and 



Il7(*)ll>0, 



d'j 
du 



are satisfied. Then for any initial estimates, any reference trajectory and any 
actuated control increments u a , the estimator errors 

||wj(i) - Wj(f)||, ie{x,y,z} 

are bounded, and the bound decreases in time for t > to- 
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Proof. Boundedness follows from the definition of the Lyapunov function V m : 

HwiW-WiWH < y/2TV m (t). 

Proposition 4.2 shows that the bound decreases in time. □ 

We do not, however, have a guarantee that the estimates will converge to 
the true values of the weights. The modified Lyapunov function V m becomes 
stationary once e m = 0, and at this point the estimation error will not be able 
to decrease further. 



4.2 Outer loop 

The task of the outer loop controller is to compute an appropriate turning angle. 
To this end we consider the outer loop nonlinear system of the form 

^ = «(cos0 t ei +sin6» t e 2 ), (20) 

where the velocity v = v(t) > is determined by the dynamics of the kite. For 
this system we obtain the following control problem: 

Problem 4.2. Find a turning angle control law such that 7 — > *y t as t — > 00. 
4.2.1 Control design 

In this section we propose a turning angle control law to solve Problem 4.2 using 
a control-Lyapunov approach. 

Recall that T t is a unit vector tangent to the target trajectory; cf Section 3.2. 
Let T := R T T t , where the operator R rotates 7 to j t along their connecting 
geodesic, so that T is an element of the tangent plane at 7. Let t := 7 x T, so 
that the vectors T and t form a basis of the tangent plane at 7. Finally, denote 

ii ■■= it - (it ■ 7)7- 

Proposition 4.3. The turning angle 6 t determined by 

cos 6 t ei + sin 9 t e 2 = T + L(t ■ 7 t L )t, (21) 
with gain L £ R + , renders the geodesic distance 

W = dist(7,7 t ) 

decreasing along the flow. 

Proof. The time-derivative of the geodesic distance is given by Lemma 6.2 in 
Appendix B: 

^ = -(7-R T 7 t )-Y} S (22) 

where Y^' is the geodesic vector pointing from 7 towards j t . Substituting (21) 
into (22), and noting that 

7t = ( 7 -T)T t , 

we obtain 



^=-Mt.Y-)(t-7^) = -^ (t - 7t±) 



2 



Note that dW/dt -> as 7^ -> 0. □ 
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4.2.2 Control performance 

Proposition 4.2.1 shows that the control law (21) renders the geodesic distance 
dist(7,7 t ) decreasing along the flow. In this sense, the control law (21) renders 
the target trajectory stable, but not asymptotically so. For the latter we would 
also need to be able to control the kite velocity v. 

5 Simulation 

In this section we report on the performance of our controller with an extension 
of a point-mass kite model commonly found in the literature. 

5.1 Model 

We present an extension of the point-mass kite model discussed in [8, 12, 13, 
15, 9]. For simplicity, our presentation uses a centripetal force to keep the kite 
on a sphere with given radius. 

We express the aerodynamic forces acting on the kite in an aerodynamic 
reference frame. Let v a = r — v w denote the apparent wind vector, where v ffi 
represents the environmental wind. Let e x be the unit vector collinear with the 
apparent wind vector v . We assume that the kite cannot roll with respect to 
the tether, and so we must have e y such that e y is parallel to r x e x . Choose 
e y to be collinear with r x e x and let e z := e x x e 9 . 

In expressing the magnitude of the forces we will make use of unitless lift 
coefficients. Let Cl ■ a i-> Cl{o) be the unitless lift curve with a :— acos(r • 
e z/ll r ll) + ^p; where 6 p is the power setting, and let Cs '■ u M- Cs(u) be the 
unitless steering force curve. We model the drag force as the sum of a zero-lift 
drag coefficient and a parabolic approximation of induced drag [2] : 

C D = C Do + kC 2 L . 

The constant k is defined by fc _1 = nAe, with aspect ratio A = b 2 /S, Os- 
wald factor e, projected surface area 5, and wingspan b. It follows that the 
aerodynamic forces are given by 

F x = - l - P \\v a \\ 2 SC D {a), 

F y = - l -ph a \\ 2 SC s {u), 

F z = -lp\\v a \\ 2 SC L (a), 

where p is the density of air. 

In our model the steering control input u causes the kite to accelerate in 
the direction of the Y-axis, which in turn causes an increase in the relevant 
component of the apparent wind vector v a . The aerodynamic reference frame 
rotates along with the apparent wind vector, and in this way the effects of 
the steering yaw moments observed in real kites are emulated. The model is 
completed by considering gravity and constraint forces to restrict the kite to a 
sphere. 
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In our simulations we choose a projected surface area of S = 11 m 2 , a wing 
span of b = 6 m, an Oswald factor of e = 0.7 and a mass of m = 1.5 kg. The 
zero-lift drag coefficient is given by Cd = 0.075. The lift curve Cl, modeled on 
data from [22] for a NACA 0018 airfoil in a flow with Reynolds number 2 x 10 6 
using B-splines, is given in Figure 6. Due to the current lack of more specific 
data, we also set Cs = Cl- This is justified by the observation that one can 
view the tips of the kite as vertically positioned airfoils, with a local attack angle 
controlled using the steering control input u. 



tS 0-0 




-90 45 90 



« [deg] 



Figure 6: Lift curve modeled on data from [22] for a NACA 0018 airfoil at Re 
= 2 x 10 6 . 



5.2 Reference trajectory 

For the tracking objective to be independent of the velocity we treat the cor- 
rected path length - as opposed to time - as the independent variable. Our 
reference turning angle is 



6[ lt ](s c ) :=A 



cos I — s c I - 1 



+ 0.834029, (23) 



where A 2.40483 is a root of the 0th Bessel function of the first kind and 
L = 4/3 is the length of the projected target trajectory. This turning angle 
generates a figure eight trajectory. It can be shown that, for the planar notion 
of turning angle, the constant A being a root of the appropriate Bessel function 
is a necessary and sufficient condition for the turning angle (23) to result in 
an L-periodic trajectory. The second constant is the initial condition which 
also corrects for the curvature of the sphere. It was found using a root-finding 
procedure. The resulting trajectory is an idealization of trajectories flown using 
the prototype Laddermill kite power plant. 



5.3 Results 

The wind vector points in the direction of the earth X-axis with magnitude 
6 m/s. The kite is tethered on a 100 m line, and the controller gains are tuned 
- using a trial and error procedure - to K = 3.0, L = 10.0, and T = 100.0. The 
steering rate is limited to 0.05 rad/s, the magnitude of the steering input to 0.12 
rad, the power setting is fixed at 9 p = 0.01 rad, and every control derivative 
estimator B-spline has 10 knots. 
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Simulation results for a run on constant line length without turbulence are 
plotted in Figures 7 and 8. Part of a trajectory with a tether reel out velocity 
equal to a third of the wind speed is shown in Figure 10. 




Figure 8: Control input for the tra- 
jectory shown in Figure 7. 



Figure 9: Mean linear tracking error 
as function of turbulence intensity a 
and correlation rate S, where <t\ = 
cia = °Vc = c Every data point is an 
average of 100 simulations. 



5.3.1 Turbulence 

In order to quantify the controller performance under the influence of turbulence 
we consider the mean linear tracking error 

1 f T 

J: =fJ o \h-7t\\dt. 

Figure 9 plots the mean linear tracking error J for a single cycle as a function of 
turbulence intensity for precise initial control derivative estimates. The turbu- 
lence is modeled according to the simplified Drydcn model [4]. The simulation 
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Figure 10: Part of a trajectory during tether reel out. 



runs at 100 Hz and every data point is an average of 100 simulations. The graph 
shows that for the investigated range of turbulence intensities the tracking error 
remains bounded by a value of 0.009. Scaling this number with the distance 
between the earth tether attachment point and the kite (100 m) we obtain a 
bound of less than one meter. 

5.3.2 Measurement noise 

Noisy measurements degrade control performance, which might conceivably lead 
to problems with adaptivity. Initial investigations suggest that measurement 
noise does not induce drift in the estimates of the control derivatives; see Fig- 
ure 11 for the evolution of a control derivative estimate with Gaussian white 
noise at the levels given in Table 1. Assessment of the exact effects of measure- 
ment noise on adaptivity performance is a topic for further investigation and is 
beyond the scope of this work. 
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Figure 11: Evolution of the 6th Y-axis control derivative estimate with Gaussian 
white noise at the levels given in Table 1. 
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Measurement 


Mean 


Standard deviation 


Position 


m 


2.5 m 


Velocity 


m/s 


0.02 m/s 


Acceleration 


m/s 2 


0.02 m/s 2 


Airspeed 


m/s 


0.9 m/s 


Attitude 


deg 


1 deg 



Table 1: Measurement noise levels. 



6 Concluding remarks 

In this paper we designed a controller for trajectory tracking control of kites. We 
projected the kite trajectory onto a sphere and used the differential-geometric 
notion of turning angle as primary tracking variable. Based on this dimension re- 
duction we derived an adaptive controller with minimal modelling requirements. 
Repeated simulations with a point-mass model show our control approach to be 
robust against turbulence. A sequel paper will detail the performance of our 
controller with multi-body simulations of a deforming kite on a flexible tether. 
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Appendix A: Geodesic curvature 

In this section we list the main results related to the concept of geodesic cur- 
vature referred to in this paper, adapted from the book by Gray [11]. We will 
not make the notions of surface and curve precise here, but only note that a 
coordinate patch is a differentiable map between an open set Wet 2 and K 3 : 
the intuition is that it associates coordinates to a part of a surface. 

In the following, J will denote a counterclockwise rotation by it/ 2 in the 
tangent plane of the surface under consideration. 

We start with the formal definition of the turning angle. Intuitively, the 
turning angle is the angle between the tangent of a curve and a reference, without 
being restricted to [— tt,tt]: 

Definition 6.1. Let M. be an oriented surface, suppose that a : (a,b) — > M. is 
a regular curve, and that X is an everywhere nonzero vector field along a. Fix 
t with a < t < b. Let 9 be a number such that 

d(to) X(t ) , . a JX(io) 

— cos — — — — + sin - 



IK*o)|| u ||X(t )|| lX(i ) 



1G 



Then there exists a unique differentiate function 9 = 0[a,X] : (a, b) — > R such 
that 8 (t ) = O n and 

M =cos^,X] M +sin^,X] M , 

at all points on the curve. We call 9[a,~X] the turning angle of a with respect 
to X determined by 9o and t . 

In order to state Liouville's theorem, which allows us to compute the turning 
angle of a curve, we need the notion of geodesic curvature. Intuitively speaking 
the geodesic curvature of a curve measures how far it is from being a geodesic, 
that is, the shortest path between two points on the surface: 

Definition 6.2. Let a : (a,b) — > M. be a curve in an oriented regular surface 
M. inM. 3 . Then the geodesic curvature of a is defined (for a <t < b) as 



a(t) ■ Ja(t) 

Hi 



Finally, we state Liouville's theorem, which relates the geodesic curvature 
of a curve to the geodesic curvatures of the coordinate lines induced by the 
coordinate patch and the turning angle of the curve. Let and (n g )2 denote 

the geodesic curvatures of the coordinate lines u M> x(u,v) and v i->- x(u,w), 
where x is a given patch. 

Theorem 6.1 (Liouville). Let Ai be an oriented surface, and suppose that 
a : (a,b) — > M. is a regular curve whose trace is contained in x(W), where 
x :U — > M. is a coordinate patch such that x„ = Jx u . Then 

Kg[a] = (Kg)l COS9 + (Kg) 2 SU1 9 + , 

where 9 = 9[ot,x u ] is the angle between a and x u . 

Appendix B: Geodesic distance on the sphere 

In this section we list the main results related to the concept of geodesic distance 
on the sphere referred to in this paper, adapted from the report by Bullo, 
Murray, and Sarti [6]. 

Definition 6.3. The geodesic distance between the points p and q on the unit 
sphere is defined as 

dist(p, q) := acosp • q. (24) 

Additionally, provided that p ^ ±q, there is a uniquely defined unit vector 
that gives the geodesic direction at the point p towards the point q: 

Definition 6.4. The geodesic vector giving the geodesic direction in a point p 
towards a point q ^ p is defined as 

Y q ._ (p x q) x p q - (q • p)p 



P ' ll(pxq)xp|| ||q-(q-p)p|r 
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Figure 12: Geodesic on the sphere. 



The time-derivative of the geodesic distance between a trajectory p = p(i) 
and a fixed point q can be expressed in terms of the geodesic direction Y£j: 

Lemma 6.1. Consider a trajectory p = p(t) on the sphere, such that p(t) never 
passes through the fixed points q or — q. Then 

|dist(p(i), g ) = -p.Y^. 

Proof. Differentiating the definition in Equation (24), we have 

^dist(p(i),q) = -^acosp-q 

p q = p (q (q-p)p) = . 

v/i-( P -q) 2 llpxqll p ' 

□ 

Let R be the rotation about pxq which maps p to q. We will now extend the 
previous lemma to cover the time-derivative of the geodesic distance between 
two trajectories p = p(t) and q = q(t). 

Lemma 6.2. Consider a trajectory p = p(t) on the sphere, such that p(t) ^ ±q(£) 
for allt > 0. Then 

|dist(p(t),q(t)) = -( P -R T q).Y^. 

Proof. By differentiating first with respect to p and then with respect to q, we 
have 

^dist(p,q) = -p.Y2-q.YP. 
Since RYJj = -YP, we can simplify the previous expression to 
d 



rft dist(p,q) = -(p-R J q).Y£. 



□ 
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